Electron-lattice coupling, orbital stability and the phase diagram of Ca2_xSr2:Ru04 
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Hartree-Fock calculations are presented of a theoretical model describing the Sr/CaRu04 fam- 
ily of compounds. Both commensurate and incommensurate magnetic states are considered, along 
with orbital ordering and the effect of lattice distortions. For reasonable parameter values, interac- 
tions disfavor orbital disproportionation. A coherent description of the observed phase diagram is 
obtained. 



The two dimensional ruthenate family Ca2_2;Sr2;Ru04 
presents a rich phase diagram Q , evolving as x is varied 
from the Fermi-liquid-triplet superconductor Sr2Ru04 to 
the antiferromagnetic (Mott) insulator Ca2Ru04. The 
origin and indeed the nature of some of the phases re- 
mains the subject of controversy. Further interest has 
been added by surface studies which indicate that 
(contrary to intuition) the Mott state is less stable on the 
surface than in the bulk. Bulk Cai.9Sro.iRu04 exhibits 
a metal to insulator transition as the temperature is de- 
creased below 150 K whereas surface sensitive probes 
place the transition at Ril25 K This unexpected re- 
sult challenges the conventional notion that correlation 
effects are stronger at surface, and requires explanation. 

In this paper, we present calculations which shed new 
light on Ca2-aiSr3Tlu04. While a prior theoretical lit- 
erature exists 0, M IE S Q I O'^'^ analysis involves new 
features including study of heretofore unexplored param- 
eter ranges, of incommensurate magnetic phases, and of 
coupling to the apical oxygen. We have also uncovered 
a previously unnoticed phenomenon: for wide and physi- 
cally reasonable parameter regimes, interactions stabilize 
the system against orbital disproportionation. 

We study a model Hamiltonian derived from a tight- 
binding approximation to the calculated band structure 
and supplemented by electron-electron and electron- 
lattice interactions: H = Hband + H^-e + He-iatt + Hiatt- 
The near-Fermi-level states are derived from Ru t2g sym- 
metry d-states (admixed with oxygen) , well described by 
a tight-binding dispersion Hjyand — H 
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Eq^ — 2f(cos kx + cos ky) — Ati cos k^ cos ky, 
= £0^'''^ - 2t2Coskx - 2t3Cosky, e'^. = 
—4^4 sin fca; sin A: J, and e\^{kx,ky) = e%^{ky,kx)- Quan- 
tum oscillation measurements imply as Sq^ = 
0.62, Eq^'^^ = 0.32, t = 0,-42,ii = 0.17, i2 = 0.30, = 
0.03 and 14 = 0.04 (eV) 9]. These band parameters give 
nxy — 1-31 and n^z — nyz — 1.35, almost 4/3. A crucial 
question, discussed in more detail below, is how the or- 
bital occupancies change as parameters are varied. The 
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We have written in Eq. ^ in the conventional notation 
in which U, U' , J, are the intraorbital Coulomb, in- 
terorbital Coulomb, interorbital exchange, interactions, 
respectively and we have omitted a pair-transfer (J') in- 
teraction which does not affect our results. For c?-orbitals, 
U, U' , and J are functions of only two of independent 
linear combinations of the fundamental atomic physics 
(Slater) parameters Fq, F2, F^. One combination involves 
the multiplet averaged interaction Fq, along with a small 
admixture of F2 and F^ and is expected to be strongly 
renormalized by solid state effects 11]; the other involves 
the combination of F2 and F4 which determines J and 
is expected to be less affected by solid state effects. The 
precise linear combinations depend on the strength of 
the ligand field and are not important here, but in cubic 
symmetry U = U' + 2 J. (In a tetragonal environment 
small corrections occur; neglected here for simplicity.) 
We therefore set U = U' + 2J and explore a range of 
U and J values, which is equivalent to changing Fq and 
bandwidth. Finally, the electron-lattice coupling and lat- 
tice distortion energy are given by 
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and Q represents a nor- 



mal coordinate of RuOg distortion (apical oxygen dis- 
placement). A and K are the electron-lattice coupling 
constant and the spring constant, respectively. 

We use the Hartree-Fock approximation to determine 
the ground state phase diagram by comparing energies 
of different phases including paramagnetic (PM), a fer- 
romagnetic (FM) state, commensurate antiferromagnetic 
states with qc = (tt, 0) (C-AFM) and qg = (tt, tt) (G- 
AFM), and an incommensurate magnetic (ICM) state 
with qic = (27r/3, 27r/3). qic is close to the momen- 
tum q* ~ (0.697r, 0.697r) where the susceptibility x has 
a local maximum due to the near nesting of quasi-one- 
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dimensional {xz,yz}-ha.TLids. Strong magnetic scatter- 
ing peaked near qjc has been experimentally observed 
in Sr2Ru04 We also considered the orbital order- 

ing (00), i. e., changing the relative occupancy of xy, 
xz and yz states. Uniform states were found, but two- 
sublattice states such as to those reported in Ref. were 
never found to minimize the energy for the parameters 
studied. 

Before presenting our results we discuss a simple one- 
atom calculation. In conventional atomic physics, one 
minimizes the energy using a state which corresponds to 
a definite electron occupancy, and finds (for n — 4) that 
the ground state corresponds to maximal spin and or- 
bital angular momentum, with excited state energies de- 
termined by the higher Slater parameters, F2, F4 i. e., by 
J. For a solid state environment, it is more appropriate to 
minimize using a density matrix corresponding to the de- 
sired mean charge density per orbital Ua and spin density 
per orbital vrio,. In this case the Fq (i. e., J7) parameter 
also contributes to energy differences. It is convenient to 
express the densities in terms of no = n^y + n-j-z + "-^2, 
Ts [defined below Eq. |Q] and T2 = ^{uxz — nyz) and 
similarly for m. For paramagnetic states (ma = 0) and 
fully polarized states (mo = 2), we find 

<r K = 4) = -i({/ - 5J)irl + rl), (5) 

Ef:'r"''in, = 2,mo = 2) = -^-{U iJ){ml + ml).{(S) 

Thus, in a solid state environment, interactions may ei- 
ther promote or inhibit orbital disproportionations; in 
the present Hartree-Fock approximation for paramag- 
netic (fully polarized) states disproportionation is sup- 
pressed for J >U/5{J > C//3). 

Figure^ shows the numerically obtained ground-state 
phase diagram without electron-phonon coupling as a 
function of U/t and J/U. At large U, two phases are 
observed; a ferromagnetic metal with a small degree of 
orbital disproportionation (FM) and an insulating an- 
tiferromagnetic phase (G-AFMOO). In the G-AFMOO 
state, the xy-orbital is occupied by two electrons and the 
other two electrons sit on the {xz,yz}-OThitals forming 
a half-filled band gapped by the AF order. Even in the 
absence of long range order, the strong correlations and 
commensurate filling of the {xz, yzj-bands would lead to 
Mott insulating behavior. This is the state proposed by 
Nakatsuji and Maeno jjj for Ca2Ru04. Because it in- 
volves orbital disproportionation the G-AF phase is only 
stable for small J/U. Turning now to the small U regime, 
we find that as U/t is increased the first phase transition 
is of second order, and is to an ICM state characterized 
by a wave vector q* ~ (0.697r, 0.697r) determined by near 
nesting of the {xz, yz} bands, x ^l^o has a strong peak 
at a much smaller 52 ^ (0.5,0.5), arising from the xy 
band, but in the present calculation the q* instability is 
slightly favored. We suggest that the IC phase is the one 
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FIG. 1: Phase diagram without electron-phonon coupling for 
band parameters corresponding to n^^y = n^z ~ Uyz =4/3 at 
U — 0. PM: paramagnetic state, FM: ferromagnetic state, C- 
AFM: C-type antiferromagnetic state [qc = (7r,0)], G-AFM: 
G-type antiferromagnetic state [qc = (7r,7r)], ICM: incom- 
mensurate magnetic state [qic ~ (3^1 3^)]- OO: orbital dis- 
proportionation with Uxy 7^ rixz.yz- Only the G-AFM state is 
insulating. Symbols connected by solid lines indicate phase 
boundaries (first order, except for the ICM/PM boundary) 
between the indicated states. Broken line: second-order phase 
transition from PM to ICM with q* determined by the suscep- 
tibility maximum. As U/t 00 the G-AFMOO-FM phase 
boundary asymptotes to J/U = 1/3 as discussed in the text. 



observed in the actual materials for 0.2 < a; < 0.5. As 
U is further increased a first order transition occurs to a 
phase which is FM or G-AFM depending on J/U. Note 
that in computing the energy of the IGM phase we ap- 
proximated the ICM vector by qjc = (27r/3, 27r/3). The 
errors due to this approximation may be estimated from 
the difference between broken and solid lines in Fig. 

We now interpret experiment in terms of our results, 
beginning with Sr2Ru04. Quantum oscillation measure- 
ments (Table 6 in Ref. j^]) directly determine both spin 
polarization and quasiparticle mass enhancement and re- 
veal a Stoner factor (enhancement of x above that due 
to the quasiparticle density of states) of 4 — 5. Within 
our approximation this corresponds to 2 < U/t < 2.3 
(weakly dependent on J/U), i. e., to a U/t « 0.8 of the 
critical value for a ferromagnetic instability (preempted 
here by the IGM state). We also find that the three di- 
mensional model has a second order PM-FM transition at 
U + 2J ^ 2.5t and for [/ -I- 2 J ~ 3t, m ~ Ifis consistent 
with observation Thus we argue that a moderate 

interaction, e. g., U ~2t and J ~ 0.6t, provides a good 
description of the Sr-ruthenates. 

Our estimate of the interaction strength implies that 
a ~ 10 % decrease in t (a change of the magnitude ex- 
pected when Sr is replaced by Ca) will drive a transition 
into an ordered state. Over a wide range of the phase 
diagram, our calculation indicates that the first insta- 
bility is to an IGM state driven by near Fermi surface 
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nesting. Indeed we find that this instability occurs at 
U /t < 2, i. e., where the Stoner factor is less than 4. 
However, beyond-Hartree-Fock-corrections are expected 
to shift the ICM phase boundary substantially, both be- 
cause we are dealing with low dimensional magnetism, on 
which quantum fluctuations have substantial effect, and 
because inelastic scattering will weaken the susceptibility 
peaks. Further, the two possible ICM wave-vectors lead 
to very similar critical J7's; beyond Hartree-Fock effects 
are expected to control whether {xz, yz}-band dominated 
g* ~ (27r/3, 27r/3) or xy-band dominated 5*2 — (0.5,0.5) 
magnetism occurs. Thus, while it is natural to expect a 
second order magnetic transition to occur as Ca is sub- 
stituted for Sr, the specifics depends on details beyond 
the scope of our calculation. We argue that the most 
natural interpretation of the data is that the transition 
observed at a: = 0.5 is the PM-ICM transition we found, 
but that the relatively larger fluctuation corrections to 
the ICM state shift the ICM critical point close to the 
FM critical point (Stoner factor at a; = 0.5 « 25). In 
the actual materials, the T = phase boundary to the 
ICM state is close to the onset of a lattice distortion [l^ 
which distinguishes the [1,1] and [1,1] directions of our 
2d lattice. In the distorted phase the two ICM vectors 
qic and qj^ = {q, —q) become inequivalent, so diagonal 
ICM order and the lattice distortion couple. A minimal 
description is 

T = TsiSjc, Sy + ^{Sjc - + ^CX\ (7) 

where Sj^fj^j^ is the magnetic order parameter associ- 
ated with qic{(fjQ), X denotes the distortion amplitude 
at 0.2 < X < 0.5 and 7 and C are a coupling constant 
and an elastic constant, respectively. Eq. shows that 
the magnetic transition will induce the lattice transition, 
consistent with the observed coincidence of the two tran- 
sitions at a; = 0.5, T = 0. Alternatively, a primary lattice 
transition could induce a first order magnetic transition. 
In either case the observed difference in transition 
temperatures follows from the strong effect of fluctua- 
tions of low-dimensional magnetism. Note, however, di- 
agonal (g, q) wave vectors are always found to be favored 
over "face" {q, 0) vectors. Bragg peaks corresponding 
to the ICM phase have not been observed. Friedt et 
al. do observe low-energy incommensurate fluctuations 
at g ~ (0.6,0) [13 but in our calculations x is always 
greater at (0.5,0.5). The coincidence of magnetic and 
structural transitions also favors diagonal order. 

Further decrease of t leads to a first order transition 
into a strongly magnetic state. With the Hartree-Fock 
approximation, the first transition is to either the C- 
AFM phase or the FM, depending on the value of J/U . 
As far as we know, neither of these phases is observed 
experimentally. To account for the direct transition into 
the G-AFM phase which is apparently observed, an ad- 
ditional coupling which favors orbital disproportionation 



must be included. A very natural possibility is coupling 
to a "Qs-type" lattice distortion (volume-preserving in- 
crease in Ru-apical-oxygen distance and decrease in Ru- 
planer-oxygen distance). Fig. El shows results obtained 
using Eq. for J/U — 0.25. We see that the range 
of G-AFM phase is increased, and for strong enough 
electron-phonon coupling, a direct transition from ICM 
to G-AFM is possible. Also within this model, a lattice 
distortion and a hardening of the phonon frequency occur 
at the transition to the G-AFM phase, again in at least 
qualitative agreement with the experiment Pressure 
experiments on Ca2Ru04 observe a transformation to a 
metallic phase with a small FM moment. Pressure in- 
creases t and K, corresponding to diagonal motion in 
Fig. [3 The FM phase shown is highly polarized in con- 
trast to that found in Ref. ^ weakly polarized FM 
phase exists as a local free energy minimum in the region 
labeled ICM; it is conceivable that pressure suppresses 
the ICM phase, revealing the FM one but this is not 
found in our theory. 

The results presented so far are obtained for a band 
structure corresponding to rixy = rixz = nyz = 4/3 at 
U — Q. Fig. O shows the evolution of orbital occupancies 
in the PM phase for different J/U and band structure 
such that at t/ = rixy = 1.4, nxz,yz — 1-3. We see that, 
as expected from the simple estimate, the orbital dispro- 
portionation is strongly enhanced by correlation effects 
if J is less than critical value Jc {— 0.2U in the Hartree- 
Fock approximation used here) and is suppressed if J is 
greater than the critical value. The strong coupling be- 
tween different types of order and the orbital occupancy 
means that for J/U < 0.2 the phase diagram is very sen- 
sitive to the xy/{xz, yz} level splitting, but for J/U > 0.2 
is much less so. 
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FIG. 2: Phase diagram as a function of U/t and X^/Kt for 
J/U = 0.25 and n^y = nxz,yz = 4/3. Solid lines: first-order 
phase transitions. Dotted line: second-order transition from 
PM to IGM. Broken line: transition to quadrupole ordering 
with rs > 0; weakly first-order because cubic terms are al- 
lowed in the free energy. 
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FIG. 3: Orbital occupancies in PM phase as functions of U jt 
for a band structure rixy = 1.4 and Uxz = riyz = 1.3 at ?7 = 0. 

Our results suggest a possible interpretation of the puz- 
zling surface experiments by Moore ei al. who find that 
the surface of Cai.9Sro.iRu04 remains metallic down to 
a lower temperature (125K) than does the bulk mate- 
rial (150K) and that both on the surface and in bulk the 
metal-insulator transition is accompanied by a Qa-type 
distortion of magnitude ~ 0.02. The surface phonon fre- 
quency of the apical oxygen mode is 10 % larger than 
bulk, implying a 20 % increase in K and thus a 20 % 
decrease in }? jKt. As seen from Fig.|21 a change of this 
magnitude makes the insulating phase substantially more 
difficult to access. A full theory requires a study of the T 
dependence, for which Hartree-Fock is unreliable. Other 
issues also remain unresolved, including the roles played 
by the frozen-in Q3 and Qo distortions occurring at the 
surface and the additional lattice strains caused by the 
bulk transition. Further investigation using dynamical- 
mean-field methods and a detailed treatment of the lat- 
tice effects would be desirable. 

Our approach is most similar to that of Nomura and 
Yamada , who performed Hartree-Fock studies of the 
model studied here, but with band parameters corre- 
sponding to slightly broader Id bands. They assumed 
U' jU = 0.5, JjU = 0.25, ('orbitally stable'), consid- 
ered only G-AFM and FM ordered phases and did not 
include electron-lattice coupling, but did investigate ef- 
fects of changes in relative level energies. In a closely 
related work. Fang and Terakura used LSDA band the- 
ory methods to investigate FM and G-AFM states Q. 
Their results, in particular their prediction of a direct 
PM-FM transition as bandwidth is reduced, correspond 
reasonably well to our results with J jV w 0.3, the ICM 
phase omitted. Very recently, the calculation was refined 
to include a LDA-I-C/ treatment of Ca2Ru04 and the ex- 
perimental structure, and results for optical absorption 
were presented. Hotta and Dagotto studied the model 
via a combination of mean-field results and numerical 
studies of small systems and showed the importance of 
coupling to the shape changes of the RuOg octahedron. 
However, their phase diagrams feature phases with com- 



plicated spatial structures which we are unable to stabi- 
lize. It is possible that the spatial structures found in 
Ref. are due to boundary effects in small size systems 
accessible numerically. Anisimov and co-workers '(3| used 
"dynamical-mean-field" and LDA-|-?7 method to study 
the U' = J = model. At strong correlation, they found 
an insulating "(2,2)" phase which is essentially the same 
as the G-AFMOO phase we find. At slightly lesser cor- 
relation strengths, Anisimov et al. find an interesting 
"(3,1)" phase in which 3 electrons are in the {xz,yz} 
bands in a Mott insulating state, while the xy band re- 
mains metallic. This work remains controversial!^ H^. 
The Hartree-Fock analogue (AF state with n^y = 1) is 
not found here: at J/U > 0.2 it is presumably disfa- 
vored by the orbital stability effects we have discussed; 
for smaller J it appears to be preempted by the C-AFM 
and ICM phases (not considered in Refs. |ESE|)- Very 
recent ARPES experiments indicate Uxz = nyz = 4/3 
[l^ . inconsistent with the (3,1) phase, but consistent 
with our orbital stability arguments of J/f/ > 0.2. 

To summarize: we have shown that moderate interac- 
tions, a reasonable J/U and electron-lattice interactions 
can account for the Ca2_2;Srj;Ru04 phase diagram, and 
argued. We argued that the materials are in the regime 
in which disproportionation is suppressed by interactions. 
Determining the susceptibility to orbital disproportiona- 
tion in other compounds and other multiplet structures 
is an important open issue. Investigation of ICM/lattice- 
distortion instability {x ^ 0.5) via Landau theory meth- 
ods and study, using dynamical-mean- field methods, of 
the G-AFM boundary are important directions. 
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